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Recent experimental and theoretical studies indicate that intramolecular energy redistribution 
(IVR) is nonstatistical on intermediate timescales even in fairly large molecules. Therefore, it is 
interesting to revisit the the old topic of IVR versus quantum control and one expects that a classical- 
quantum perspective is appropriate to gain valuable insights into the issue. However, understanding 
classical phase space transport in driven systems is a prerequisite for such a correspondence based 
approach and is a challenging task for systems with more then two degrees of freedom. In this work 
we undertake a detailed study of the classical dynamics of a minimal model system - two kinetically 
coupled coupled Morse oscillators in the presence of a monochromatic laser field. Using the technique 
of wavelet transforms a representation of the high dimensional phase space, the resonance network 
or Arnold web, is constructed and analysed. The key structures in phase space which regulate the 
dissociation dynamics are identified. Furthermore, we show that the web is nonuniform with the 
classical dynamics exhibiting extensive stickiness, resulting in anomalous transport. Our work also 
shows that pairwise irrational barriers might be crucial even in higher dimensional systems. 



I. INTRODUCTION 

The influence of intramolecular vibrational energy re- 
distribution (IVR) on the control of quantum dynamics 
has been a subject of some debate ever since the first 
lasers gave hope to the dream of doing mode-specific 
chemistry [TH5]. In fact, arguing that the energy flow 
process is typically faster than the reaction timescales [6] 
and essentially statistical, researchers identified IVR as 
the main culprit which spoils the idea of mode-specific 
control[7]. Consequently, most of the quantum control 
methods shifted their focus away from explicitly control- 
ling or beating the IVR process. Instead, the various 
control strategies [SHIP] concentrated on manipulating 
the quantum interferences, for which a prior knowledge of 
the field-free IVR mechanism is supposedly not required. 

The observations above, however, need to be reap- 
praised in light of the recent developments in our 
understanding of IVR. Both experiment al[lTHl3] and 
theoretical 1 14ffTo] studies indicate that IVR at interme- 
diate times is nonstatistical even in moderately sized 
molecules. Specifically, from a molecular state-space 
perspective [T7HT9] . IVR is anisotropic with local, as 
opposed to global, density of states playing a crucial 
role and significant deviations from the Fermi golden 
rule regime in the form of power law decays are to be 
expected [201I2T]. These studies certainly give rise to new 
hopes for controlling the IVR process itself and there has 
been some promising work [2"2Tt24| in this direction. How- 
ever, in light of the recent advances, there is also a case to 
be made for reopening the old debate - can mechanistic 
insights into the energy flow dynamics lead to better con- 
trol strategies? It is not entirely unreasonable to expect 
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that a subtle interplay between the IVR dynamics and 
the controlling field must occur in many of the quantum 
control studies j25j [26]. Hence, is it possible to envisage 
situations wherein IVR might actually aid, rather than 
thwart, the control process? An interesting "inverse" 
question can be asked at this stage. What features of 
a control field, as obtained through optimal[27] control 
theory, are a direct consequence of the IVR dynamics? 
There is a striking parallel here to the more contempo- 
rary issue of the search for decoherence-free subspaces 
in quantum control of open systems [28 . A mechanis- 
tic understanding of the interplay between decoherence 
and control of a system is proving to be invaluable in a 
wide range of areas. Therefore, although answering the 
inverse question posed above is far from easy, one hopes 
that a careful study of model systems can lead to a finer 
understanding of the role of IVR in quantum control of 
molecular processes. 

In order to uncover the mechanism and influence of 
IVR on quantum control one can adopt a purely quan- 
tum, purely classical, or a semiclassical approach. How- 
ever, there is little doubt that detailed insights come from 
studying and comparing all of the approaches. In particu- 
lar, although quantum studies of fairly high dimensional 
systems are now possible, a key reason for taking the 
classical-quantum correspondence perspective has to do 
with the fact that classical dynamics provides the correct 
baseline to gauge the extent to which quantum effects 
like interference and tunneling are important for con- 
trol. Such a classical-quantum correspondence perspec- 
tive has proved useful not only in the context of IVR and 
control but also in many other areas like driven atomic 
systems [21], mesoscopic systems 30J and trapped cold 
atoms [31] . Various studies on systems with two or fewer 
degrees of freedom (N < 2) have led to the identifica- 
tion of key phase space structures that regulate IVR and, 
more importantly, their quantum fingerprints. Thus, 
the role of chaos, Kolmogorov-Arnold-Moser (KAM) bar- 
riers, partial barriers like broken separatrices|32j and 
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cantori 33 , and phase space bifurcations [33] have been 
extensively studied [35| [36]. 

Nevertheless, there are two major obstacles in extend- 
ing the correspondence approach to systems with more 
than two degrees of freedom. The first obstacle is the dif- 
ficulty in visualizing the global dynamics in high dimen- 
sional phase spaces. The second obstacle has to do with 
the fact that the nature of phase space transport changes 
dramatically in going from two to three degrees of free- 
dom. For example, KAM tori are no longer true barriers 
in the phase space and there are no obvious generaliza- 
tions of the notion of cantori in N > 3. In addition, 
Hamiltonian systems typically exhibit "stickiness" [30] 
i.e., classically chaotic trajectories spending significant 
amount of time near regular structures or their remnants 
in the classical phase space. Stickiness leads to long time 
dynamical correlations and is key to understanding the 
nature of phase space transport. Important insights into 
the universality of stickiness in Hamiltonian systems are 
starting to emerge [4TH17] but their quantum manifesta- 
tions, if any, are mostly unexplored. These obstacles, 
technical as well as conceptual, have been largely re- 
sponsible for a waning of interest among researchers in 
studying and utilizing the insights afforded by a classical- 
quantum correspondence approach. 

In the last few years, however, there has been a rapid 
progress in our ability to investigate phase space trans- 
port in such high dimensional systems. Partly, this has 
to do with the recent exciting work |48ff50] in transition 
state theory - construction of locally recrossing free di- 
viding surfaces in phase space and dynamics near sad- 
dles on potential surfaces for N > 3 systems. In turn, 
such studies have led to researchers extending the con- 
cept of a transition state for driven systems [5T] with ob- 
vious relevance to quantum control. On the other hand, 
attempts are also being made[52-54 to connect the na- 
ture of global phase space transport with the structure 
of the nonlinear resonance networks, also known as the 
Arnold web. In particular, the Arnold web is central for 
a classical description of the IVR dynamics. Although 
a precise correspondence between the classical dynamics 
on the Arnold web and the quantum dynamics in the 
molecular state space is lacking, there is sufficient rea- 
son to expect that the classical insights will prove cru- 
cial in unearthing the mechanisms of IVR in a molecule. 
Will the classical insights also help in controlling the 
intramolecular energy flow itself or quantum control of 
processes influenced by IVR? At present this question 
remains unanswered. These recent advances combined 
with some of the early pioneering work collectively moti- 
vate the present study, which is a first step in answering 
the question posed above. 

In this work we take a minimal model of two kinetically 
coupled Morse oscillators, describing two local modes, 
driven by a monochromatic field to understand the in- 
terplay of IVR and control in molecular dissociation dy- 
namics. More specifically, we consider parameters appro- 
priate for the HCN (hydrogen cyanide) molecule. HCN 



has been a benchmark "small" molecular system for sev- 
eral quantum control studies [55 -58 whose primary, and 
rather challenging, aim is to dissociate the stronger C=N 
bond while keeping the weaker C—H mode "quiet" . How- 
ever, in the current study we do not attempt the selective 
breaking of the C=N bond and instead study the C—H 
bond dissociation. The focus here, in this N = 2.5 de- 
grees of freedom system, is on identifying the main classi- 
cal phase space structures that regulate the dissociation 
dynamics. In Sec. [H] we give a brief description of the 
model Hamiltonian. The details of the construction of 
the Arnold web are discussed in Sec. |III| and the impor- 
tant features in phase space arc identified. Interestingly, 
it is observed that the crucial regions in the web corre- 
spond to classical phase space structures with highly irra- 
tional frequency ratios. These results are reminiscent of 
the early work by Martens, Davis, and Ezra[59 wherein 
the importance of pairwise noble barriers to IVR in the 
OCS molecule (N 
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3) was highlighted. In Sec 
study the role of two important regions ("hubs") in the 
resonance network. One of them, the "dissociation hub" , 
acts as a gateway for dissociating trajectories and pos- 
sibly of some interest in the context of transition state 
theory. In the other hub, the "noble hub", we establish 
the extensive stickiness of the dynamics with significant 
implications for transport in the phase space. Section [TV] 
concludes. 



II. PRELIMINARIES: MODEL HAMILTONIAN 
AND MOTIVATION 

We consider two kinetically coupled Morse oscilla- 
tors being driven by a monochromatic field. The corre- 
sponding classical Hamiltonian, considered in an earlier 
work [56] as well, is given by 
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where the matter part of the Hamiltonian is 
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and the matter-radiation interaction in the dipole limit 
is taken to be 
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The form of the dipole function^^ is given by the ex- 
pression 



L ( x ) = e ^ x+x ^J2Mx + x e y 
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and the values of the various parameters can be found in 
the earlier work[56. 

As indicated In Eqn. [T] the various parameters of H 
are appropriate for the HCN molecule. Thus, Hqh and 
Hcn are Morse oscillator Hamiltonians for the anhar- 
monic vibrational displacements of the C— H (x) and 
C=N (y) bonds respectively with x e being the equilib- 
rium bond length of the CH mode. The reduced masses 
of the CH and the CN modes are denoted by Men and 
hlcN respectively and Mc denotes the mass of the car- 
bon atom. Note that the number of bound states of 
the CH and the CN oscillators are about Nb = 27 and 
Nb = 80 respectively for the given parameters. The vi- 
brational local modes are directly coupled to each other 
via the kinetic coupling Tcoup and is responsible for IVR 
in the system. In addition, an indirect coupling is also 
present between the two modes due to the interaction 
with the laser field via Hj. 

Throughout this work the various parameters corre- 
sponding to the Morse oscillators and the driving field are 
all in atomic units. The driving field strength A^ = 0.009 
(~ 3 TW/cm 2 ) and frequency uj f = 0.011 (~ 0.5 fs" 1 ) 
are chosen to avoid any significant excitation of CN bond, 
since we are interested in studying the effect of driving 
field only on the CH mode. 

A few comments on the choice of the model Hamilto- 
nian and the form of the matter-field coupling is appro- 
priate at this stage. Firstly, the form of Hq has a long 
history in chemical physics starting with the beautiful 
work by Thiclc and WilsonJSD] and the extensive studies 
by Oxtoby and Ricejl.Bi^ on the connection between non- 
linear resonance overlaps and IVR. This choice is simple 
enough to use the analytically known action-angle vari- 
ables for a Morse oscillator but has the basic features of 
interest to us - IVR, bifurcations, and a mixed regular- 
chaotic phase space. Secondly, as mentioned in the intro- 
duction, our goal here is to mechanistically understand 
the role of IVR in regulating the dissociation dynamics 
and identify the key phase space structures. In this initial 
study, therefore, we choose to avoid complications from 
other competing processes that may hinder a clear inter- 
pretation. Thus, we choose a simple monochromatic form 
for the driving field and it is important to note that this 
model already has a five-dimensional phase space. For 
the same reason we choose to couple only the CH mode 
to the field. For the weak driving fields (the ionization 
threshold for HCN is around 100 TW/cm 2 ) and the spe- 
cific initial states considered here (see next section) this is 
reasonable. Moreover, even upon coupling the CN mode 
to the field, our study indicates no significant change in 
the qualitative nature of the phase space structures dis- 
cussed later. This minimal model, therefore, has the ba- 
sic ingredients i.e., IVR , a field attempting to dissociate 
a bond, and an essential degree of dynamical complex- 
ity. In addition, such a model represents the next level 
of difficulty in extending the insights gained from studies 
on one dimensional driven Morse oscillator system [55] to 
systems with higher degrees of freedom. Note that the 



parameters are representative of HCN but the model is 
still a reduced dimensional one since the bending mode 
is not taken into account. It is known from a recent 
work [55] that the bend mode is important to the dissoci- 
ation dynamics. However, understanding the phase space 
transport in the current reduced dimensional model is a 
prerequisite for attempting the more challenging case of 
including all three modes of the molecule. 



A. Motivation 

In order to motivate the current study and connect 
to the discussion in the introduction we start by in- 
vestigating the dissociation dynamics of the zeroth or- 
der state \iich ,ncN) — |8, 3). This specific choice of 
state has to do with the fact that in the absence of laser 
field the classical phase space exhibits several interesting 
structures at the corresponding zeroth-order energy. In- 
deed, the Poincare surface of section shown in Fig. [ija) 
exhibits mixed regular-chaotic regions with prominent 
f!cH : ^cn = 1:1 and flcn '■ ^cn — 3:2 nonlinear res- 
onances interspersed in the chaotic sea. In addition, the 
Husimi distribution [64 of |8, 3) superimposed on the sur- 
face of section in Fig. Qa) shows that the state is mostly 
localized in the VLqh ■ ^cn = IT resonance region in 
the phase space. Note that |8,3) is not an eigenstate of 
the matter Hamiltonian Hq and thus the Husimi distri- 
bution does spread across the border between the main 
1:1 island and the chaotic sea. For this typical N = 2 
KAM system, the regular-chaos border has a hierarchical 
structure which is expected to result in sticky dynamics. 
In particular, a very high order resonance island is visible 
in the border region in Fig. [lja) which suggests the exis- 
tence of cantori in the phase space. Such barriers, partial 
or otherwise, limit the extent of classical IVR in the sys- 
tem. Specifically, in the undriven subsystem (Xp = 0), an 
analysis along the lines of the earlier work[61 by Oxtoby 
and Rice reveals that the overlap of low order resonance 
zones is not extensive enough at the energy of interest 
to result in dissociation. Thus, IVR alone is incapable of 
dissociating the CH bond. This is confirmed in Fig. [ljb) 
where the classical dissociation probability is zero for the 
entire time duration of interest. We stress at this junc- 
ture that any other initial state is equally a good choice 
as long as the generic criteria of a mixed phase space 
and partially localized Husimi distribution are satisfied. 
For instance other states such as |0, 15) and |4, 8) yield 
similar qualitative insights as the one chosen here. 

We now introduce the driving field with parameters 
mentioned above and note that ujp = 0.011 is chosen 
so as to be off-resonant with respect to the zeroth order 
state of interest. To begin, we decouple the CH and CN 
modes by turning off the kinetic coupling term T coup in 



Eq. 2c and drive the CH mode. Classical trajectories are 
propagated by integrating the Hamilton's equations of 
motion up to a final time 50001V with Tp = 2?r /ujp ~ 14 
fs being the field period. Trajectories are considered [63] 
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FIG. 1. (a) Husimi distribution of the initial state |8, 3) superimposed on the corresponding CH-mode Poincare surface of 
section in the absence of the driving field, (b) The classical dissociation probability P% for the CH bond in the presence of 
the driving field (triangles) with strength Xf = 0.009 and frequency lof = 0.011. Note that the CH bond does not dissociate 
until 50007V in the absence of the driving field (squares, IVR alone) or when the kinetic coupling between CH and CN mode 
is switched off (circles, field alone). 



to be dissociated if the displacement of the bond of in- 
terest exceeds a certain threshold, set here to be 15 a.u. 
for the CH bond dissociation. For the chosen driving 
field strength it is clear from Fig. [ljb) that again there 
is negligible dissociation of the CH bond [66 . Hence, the 
mode-mode coupling i.e., IVR in the absence of the field 
or the field in the absence of IVR are incapable of dis- 
sociating the CH bond. However, and this is the key 
observation, the weak driving field with nonzero mode- 
mode coupling does lead to large dissociation probability, 
as seen in Fig. [TJb) (triangles) and poses an interesting 
question. What is the precise mechanism, undoubtedly 
involving a subtle interplay between IVR and the driving 
field, by which the initial state |8, 3) undergoes dissoci- 
ation? In the following section we construct the Arnold 
web|67j and attempt to answer the question from the 
perspective of transport in the classical phase space. 



III. TIME-FREQUENCY ANALYSIS AND 
VISUALIZING THE PHASE SPACE 

The example summarized in Fig. [I] motivates us to 
understand the phase space transport in the system 
of interest. However, the phase space in this case 
is five dimensional and not easily visualized. On the 
other hand, recent studies on IVR from time dependent 
perspective [53[ [551 have established the importance of the 
resonance network or Arnold web i.e., the various reso- 
nances and their disposition in the phase space. There- 
fore, following the dynamics of an initial state on the 
resonance network should reveal the dominant frequency 
lockings experienced by the state as a function of time. 



Consequently, one can obtain valuable information re- 
garding the local structures in the phase space which reg- 
ulate the dissociation dynamics in the presence of laser 
field. Our interest here is to test if such an approach 
would indeed provide the mechanistic insights. 

To begin with one needs to construct the Arnold web 
and then interpret the dynamics of the system. The 
implementation aspect has been dealt with by several 
studies over the last decade or so. Starting with the 
early pioneering work by Martens, Davis, and Ezra|59j 
on the use of windowed Fourier-transforms and a lit- 
tle later by Laskar[6S;, several methods like fast Lya- 
punov indicators [37) , mean exponential growth of nearby 
orbits 38], and frequency modulation indicator [39 have 
been proposed. However, the interpretation of the classi- 
cal dynamics is not so straightforward for nonintegrable 
systems for reasons mentioned in the introduction. In the 
next subsection we utilize a wavelet based technique, pro- 
posed by Arevalo and Wiggins [55], to construct the res- 
onance network for the system of interest. Subsequently, 
we attempt to provide an interpretation of the dynamics 
based on the network and highlight the role of IVR in 
presence of an external driving field. 



A. Constructing the resonance network 

Since the exact action-angle variables (J, ff) for a Morse 
oscillator are known one can explicitly write down the 
zeroth-order nonlinear frequencies. Thus, for the CH os- 
cillator one has 
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with u>ch being the harmonic frequency of the CH mode. 
A similar expression for the CN oscillator can be writ- 
ten down. The general resonance condition involving the 
two oscillators and the driving field frequency can be ex- 
pressed as 

fcifi&L + fcafi&k + k s u; F = (6) 

with (ki,k2,ks) being a set of mutually prime integers. 
The order of a given resonance is Y) j \ kj | and the set 
of all resonances up to a given order forms the Arnold 
web at that order. Figure [2] shows the web of resonances 
up to a total order of five. Pure mode-mode (binary, 
k a = 0) resonances appear on the web as lines with pos- 
itive slopes and do not involve the field. In contrast, 
ternary resonances (fci, ki > and k^ < 0) involve both 
the modes and the field and appear as lines with negative 
slopes. Intersection of the binary and ternary resonances 
result in mode-field resonances (vertical and horizontal 
lines) and such intersections form resonance junctions. 
Some of the resonance junctions, for example at (1,1), 
(1/2,1), and (3/2,1/2) in Fig. [2] are clearly visible. The 
resonance junctions are expected to be important for the 
phase space transport. However, such a "static" con- 
struction cannot a priori highlight the dynamically sig- 
nificant regions of the web except in the limit of small 
perturbations [551 155"] . The system of interest here is cer- 
tainly far from near-integrability (cf. Fig.[lja)) and hence 
the static web can only act as a guide. In particular, 
irregular trajectories can have several frequency compo- 
nents and even the dominant among them evolve non 
trivially with time. It is hence necessary to construct the 
"dynamical" web from the actual nonlinear frequencies 
(£Ich (t), QcN(t)) determined as a function of time. 



Classical trajectories are computed over a time inter- 
val [0, T] starting with the zeroth-order actions corre- 
sponding to the initial state of interest. The dynamical 
phase space function zi{t) = xi(t) + ipi(t) correspond- 
ing to the Z th mode is subjected to a continuous wavelet 
transform 69 

/°° ft — b\ 

Zi(tW (-— jdt, (7) 

with a > and real b. As in the previous studies [53 ] 154 ] 
155] we choose a Morlet- Grossman form for the mother 
wavelet 

yj(t) = _J_ e ^t e ~S/2*^ (8) 
ay'2-K 

with A = 1 and a = 2. The wavelet transform in Eq. [7] 
yields the frequency content of zi (t) over a time window 
around t = b. Typically, for an integrable or regular mo- 
tion there is a single dominant frequency characterizing 
the phase space structure over the entire time. How- 
ever, an irregular trajectory is associated with several 
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FIG. 2. Arnold web for the system of interest showing all 
the resonances of order ^ \kj\ < 5, satisfying the condition 
in Eq. [5] as a frequency ratio space. As a guide the lines are 
shown with thickness inversely proportional to the order of the 
resonance. Vertical and horizontal lines correspond to mode- 
field resonance whereas positively sloped lines correspond to 
mode-mode resonance. The dotted lines show some of the 
possible pairwise noble barriers. See text for discussions. 



frequency components. In principle it is important to 
take into consideration all the frequency components of 
the time evolved classical trajectories. However, here we 
work with the dominant frequency at a given time which 
is computed by determining the scale (a, inversely pro- 
portional to frequency) which maximizes the modulus of 
the wavelet transform, i.e., fli(t = b) — max a \L^zi(a, b)\. 
This yields the required nonlinear frequencies Q(t) and 
the dynamics is followed in the frequency ratio space [59] 
(^Cir/^F, ^cn /^f)- Note that previous works suggest 
that the web generated using the dominant frequencies 
is still capable of providing important information on the 
nature of the classical dynamics |53l [54]. 

The frequency ratio space is divided into cells and the 
number of times that a cell is visited by all the trajecto- 
ries is presented as a density plot. We further normalize 
the highest density to one for convenience. For the sake 
of clarity the resulting network is shown on a logarithmic 
scale in order to highlight the significant regions. 



B. Interpretation of the dynamical resonance 
network 

The resulting density plot for the state |8, 3) is shown in 
Fig. [3] in the frequency ratio space (VLch/^f^cn /^f)- 
Comparing this dynamical frequency ratio space to the 
static construction shown in Fig. [2] reveals several crucial 
aspects of the dissociation dynamics. First and foremost, 
the heterogeneous nature of the dynamical frequency ra- 
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tio space indicates that the phase space dynamics is 
nonuniform. Thus, there exist dynamically significant 
regions in the phase space which, based on the earlier 
work [52] by Shojiguchi et al. on a different model sys- 
tem, are expected to strongly influence the dissociation 
dynamics. It is worth mentioning that several regions 
of the static web are not dynamically relevant for the 
specific initial state and timescale of interest. Similar 
observations have been made [68 earlier in the context of 
coupled standard maps. Secondly, a relatively narrow re- 
gion of the frequency ratio space around £Icn/<^f ~ 0.9 
corresponds to the dissociation of the CH mode (char- 
acterized by low values of Qch/^f)- Thirdly, a careful 
look at the web in Fig. [3] indicates the presence of three 
distinct regions in the phase space and below we focus 
on the nature of these regions. 

The first prominent region in Fig. [3] corresponds to 
the £Ich ■ ^cn = 1:1 nonlinear resonance. Note 
that the Husimi distribution of |8, 3) in Fig. [ija) does 
exhibit substantial localization in the large 1:1 island. 
However, the enhanced density observed in Fig.|3]is non- 
trivial since the invariant tori associated with the 1:1 
resonance are no longer true barriers in presence of the 
driving. As seen later, the classical trajectories do leave 
the resonance zone but repeatedly get captured into the 
resonance with large residence times. Consequently, the 
enhanced density near the 1:1 resonance has serious im- 
plications for the nature of the phase space transport. 
A second important region of the frequency ratio space 
is marked as "NH" (noble hub) in Fig. [3] with most of 
the activity occurring around a hub consisting of sev- 
eral mode-mode irrational (noble) frequency ratios and 
high-order resonances. The noble numbers are highly ir- 
rational numbers of the form (n + n'jg)/(m + m'"/ g ) with 
nm' — ran' — ±1 and 7" 1 = (v5 + l)/2 being the golden 
ratio. High density around Qch — &cn resonance and 
the noble hub indicates continuous interaction between 
CH and CN modes via IVR that seems to regulate the 
dissociation dynamics of the CH bond in the presence of 
the driving laser field. An equally important region of the 
frequency ratio space, marked as "DH" (dissociation hub) 
in Fig. |3j is in the vicinity of the Qch/^f — lg — 1 irra- 
tional ratio. Interestingly, the narrow dissociation zone 
mentioned above lies mainly to the left of this DH region. 

The enhanced density at the NH and DH hubs suggests 
an intriguing possibility - regions with pairwise noble fre- 
quency ratios seem to be acting as barriers to transport 
in higher degrees of freedom systems. At this stage, it 
is worth recalling an earlier classic study by Martens, 
Davis and Ezra[5j|I on the OCS molecule wherein it was 
conjectured that pairwise noble barriers in phase space 
might play a crucial role in the dynamics of systems with 
more than two degrees of freedom. A recent detailed 
study[7D] by Paskauskas, Chandre, and Uzer clearly es- 
tablishes that two-dimensional invariant tori can act as 
phase space bottlenecks to intramolecular energy flow in 
OCS. Therefore, the observation of the DH and NH hubs 
in the current work lends further support to the notion 



that lower dimensional phase space structures may act as 
traps and barriers for the classical dynamics in systems 
with high degrees of freedom. We now turn to a closer 
analysis of the nature of these dynamically significant 
regions of the web. 
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FIG. 3. Dynamical frequency ratio space (Arnold web) corre- 
sponding to the initial state \tlch, ticn) = |8, 3), constructed 
via time-frequency analysis of an ensemble of 5000 trajecto- 
ries propagated from r = to r = 5000 2>. Note that the 
normalized density is plotted on the logarithmic scale in or- 
der to highlight the sparse yet key regions in the phase space. 
Two key regions denoted DH and NH (enclosed by squares), 
the location of the initial state (black circle), and the CH 
dissociation region are also indicated. 



C. Role of the dissociation hub 

The global phase space dynamics shown in Fig. [3] in- 
volves an ensemble of trajectories with constant actions 
but varying angles. In order to understand the role of the 
DH and NH hubs it is useful to distinguish between var- 
ious dynamical classes of trajectories. Specifically, the 
classification in terms of regular, chaotic but nondisso- 
ciating, and chaotic dissociating trajectories can be ob- 
tained from the computed nonlinear frequencies. 



In Fig. |4j[ a )j a representative chaotic trajectory with a 
very large dissociation time is shown in the frequency ra- 
tio space. It is evident from the figure that the trajectory 
jumps between the important hubs and resonances in the 
web. However, interestingly, at a specific time the trajec- 
tory crosses the DH hub and promptly dissociates. We 
observed this to be a feature that is common to several 
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FIG. 4. Examples of different classes of dissociating chaotic trajectories and the importance of DH and NH regions, (a) A 
trajectory shuttles between resonance zones and eventually dissociates at 3300 Tf- (b) A fast dissociating trajectory which 
rapidly dissociates by 300 Tf- Note that the trajectories in (a) and (b) cross the DH region, indicated in Fig. [3] and then 
dissociate, (c) Frequency ratio space for |8, 3) constructed by considering the trajectories that do not dissociate in presence of 
the laser field until the final time of 5000 Tf- Note that in (a) and (b) the points are plotted at intervals of 6.5 Tf and 2.5 Tf 
respectively for clarity. 



such trajectories which dissociate over long timescales. In 
contrast, a rapidly dissociating chaotic trajectory shown 
in Fig. [4jb), undergoes very different dynamics enroute 
to dissociation. On comparing with the dynamics in 
Fig. [4ja) , there is a marked decrease in the jumps be- 
tween the ftcH ■ Qcn = 1:1 region and the NH region. 
Nevertheless, prompt dissociation again seems to occur 
once the trajectory crosses the DH region at a specific 
time. On the other hand, a regular trajectory (not shown 
here) remains locked in the Qch ■ &cn = 1:1 resonance 
and a chaotic but non-dissociating trajectory jumps be- 
tween the flcH ■ &cn = 1:1 resonance and the noble 
hub over the entire time duration. More importantly, the 
nondissociating trajectories, regular and chaotic, do not 
cross the specific DH region in the frequency ratio space. 

Is crossing of the dissociation hub a generic feature of 
all the dissociating trajectories? In order to answer this 
question, in Fig. we show the frequency ratio space 
from a data set containing only non-dissociating trajec- 
tories. Interestingly, Fig. |4jc) shows that the density in 
the vicinity of the DH hub diminishes whereas the density 
around the NH region remains unchanged. This clearly 
indicates that the high density at the dissociation hub in 
Fig. [3] is a specific characteristic of dissociating trajec- 
tories. Therefore, the DH region is like a "gateway" for 
dissociation since the frequency ratio space is also rather 
sparse to the left of the DH region. 



D. Role of the noble hub 

It is clear from the web shown in Fig.[3]that 1:1 mode- 
mode resonance and the noble hub consisting of several 
irrational frequency ratios, and their rational approxi- 
mants in terms of high order resonances, are dynami- 
cally the most visited regions of the web. The enhanced 
density in the NH region could be due to two main rea- 



sons. Firstly, this might be because of frequent visits of 
every trajectory in the specific region of the frequency 
ratio space resulting in the observed enhanced density. 
Secondly, the trajectories could be genuinely sticky and 
upon entering the region could be trapped for a consider- 
able time. Such stickiness of individual trajectories can 
also lead to enhanced density. Establishing whether the 
NH region is sticky or not is crucial since the existence of 
sticky regions in phase space imply long time dynamical 
correlations which, in turn, have important implications 
for the transport mechanism [40 j. Note that Fig. Qc) 
already hints at the possibility of NH being a sticky re- 
gion. In addition, our calculations (not shown here) of 
the mean square displacement of the action Jch 

((AJ ch ) 2 )(t) = ([J C h(t) - Jch(0)} 2 ) ~ (9) 

with the average being done over nondissociating trajec- 
tories, show /3 < 1 and hence subdiffusive behaviour. 



Further understanding into the nature of the NH re- 
gion can be obtained by analysing the residence times of 
trajectories entering the region. The approach employed 
to obtain the residence time distribution is as follows. We 
define boxes in the high density regions of the frequency 
ratio space as shown in Fig. [3] During the dynamics 
the time at which the trajectory enters and then leaves 
this region are denoted by n and T2 respectively. The 
residence time r r = t-i — t\ for one such event of the tra- 
jectory residing in the NH region is noted. There may 
be several events of such kind for a trajectory and out of 
all such events we pick the longest r r in the region of in- 
terest. We then construct the distribution of the longest 
r r for the entire ensemble of trajectories for fixed initial 
actions (corresponding to the initial state of interest) and 
varying initial angles. 




FIG. 5. Distribution <j> of longest residence times in specific regions of the frequency ratio space corresponding to the initial 
state |8, 3). The panels from left to right show cf> for the NH region, the 1:1 mode-mode resonance region and the DH region 
respectively. 



The distribution of the longest residence times in the 
NH, CIch = &cn resonance and DH regions are shown in 
Fig. [5] Clearly, the resonant and the NH regions do sup- 
port trajectories with fairly long residence times. There- 
fore, the enhanced density in the NH region of the web 
in Fig. [4] arises due to extensive stickiness and supports 
the subdiffusive behaviour of ((AJch) 2 )(t) mentioned 
above. In comparison, the DH region is far less sticky 
and hence in agreement with it being interpreted as the 
dissociation gateway. Although not shown here, studies 
on other initial states do highlight similar regions in the 
frequency ratio space. The differences in the dissocia- 
tion dynamics, however, arise from the relative extent of 
trapping in the various regions. The results in this sec- 
tion imply that at least on the classical level it is not 
appropriate to describe the dissociation as being statis- 
tical. The three main hubs in the classical phase space 
play a crucial role in the dissociation mechanism. 

As the key regions in the phase space have been identi- 
fied, the next logical question is whether we can discover 
the mechanism of dissociation. In order to answer this 
question we also need the precise mechanisms, in terms 
of appropriate phase space structures, by which trajecto- 
ries are transported from one region to the other. In par- 
ticular, one needs to understand the transport near reso- 
nance junctions - a demanding[7T] but important task for 
the future. Nevertheless, an idea of the different mech- 
anisms can be obtained from the dynamical frequency 
ratio space. The jumps in Fig. |4|a) happen to be along 
lines with negative slopes. On comparing to the static 
web in Fig. [2] it seems that the ternary resonances are 
playing a crucial role. Although there are several high 
order mode-mode resonances in the vicinity of the NH 
region, it is useful to remember that the CH mode does 
not dissociate in the absence of the field. Clearly, signifi- 
cant excitations in the CN mode are required in order to 
cause extensive IVR and evidence for the same can be ob- 
served in Fig. [4] as visitations of the trajectory to regions 
with low values of Qcn- As a rough estimate, a value 
of Hcn/^f ~ 0.7 corresponds to the unperturbed action 
value of Jcn ~ 20. Thus, the ternary resonances along 



with the 1:1 resonance must be enhancing the IVR be- 
tween the two modes, ultimately resulting in dissociation. 
However, Fig. [4] unequivocally shows that the pathways 
of the short and long time dissociation trajectories are 
different. 



IV. CONCLUSIONS, CHALLENGES AND 
FUTURE DIRECTIONS 

In this work we have studied the classical dynamics of 
a system with N = 2.5 degrees of freedom to gain in- 
sights into the dissociation dynamics in presence of IVR. 
Our work emphasizes the utility of time-frequency anal- 
ysis in understanding the nature of phase space trans- 
port in such high dimensional phase spaces. In partic- 
ular, we have constructed and identified regions in the 
resonance web i.e., in the full phase space, which con- 
trol and regulate the classical dissociation dynamics of a 
model triatomic molecule. Three key regions in the web 
regulate the dynamics. Of these, a region called as the 
dissociation hub acts as a gateway to dissociation and 
the region called as the noble hub, corresponding to the 
overlap of several high order mode-mode nonlinear reso- 
nances and pairwise irrational frequency ratios, leads to 
extensive stickiness. The combined influence of the disso- 
ciation and noble hubs in the phase space results in the 
phase space transport being anomalous. Thus, despite 
being chaotic, certain trajectories survive for very long 
time without undergoing dissociation. Moreover, the fre- 
quency ratio space dynamics clearly shows the existence 
of different dissociation pathways for prompt and late 
dissociating class of trajectories. 

However, the present study represents only a start to- 
wards, what we believe, are important questions that 
need to be answered in order to understand and influence 
the mechanism of dissociation. We conclude by mention- 
ing a few of these issues. 

1. What is the precise mechanistic difference between 
the fast and slow dissociating trajectories? Some 
hints do emerge from the present work. However, 
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a systematic study of the role of the various res- 
onance junctions and the phase space structures 
which mediate transport between them is necessary 
to make further progress. In this regard it is im- 
perative to construct the full Arnold web and 
study the time evolution of the web. 

2. Is it possible to relate the notion of the dissociation 
hub as a gateway to the concept of transition states 
in driven systems |51|? This requires one to deter- 
mine the invariant manifolds in phase space using 
the technique of time-dependent normal forms and 
thereby distinguish between long time chaotic dis- 
sociating and nondissociating trajectories. Such a 
study would provide the much needed connection 
between the dynamics in the frequency ratio space 
and the powerful approach of transition state the- 
ory. 

3. Do the classical phase space structures and the re- 
sulting mechanisms survive quantization? This is 
an important issue and the studies reported here 
provide a firm baseline for establishing a clear 
classical-quantum correspondence perspective into 
the role of IVR in quantum control. Preliminary 
computations of the quantum dissociation proba- 
bilities suggest that changes in the phase space dy- 
namics caused by external driving fields, especially 
stickiness, do influence the quantum transport as 
well. In particular, there are indications that the 
various hubs are more restrictive in the quantum 
case. 

4. Given that the key phase space structures regulat- 
ing the dissociation dynamics have been identified, 
is it indeed possible to design control schemes by 
interfering with the specific phase space structures? 
Such a local phase space approach towards control 



has been explored recently for systems with lower 
degrees of freedom [72} [73] . It remains to be seen 
if the insights gained from the current work can 
be utilized for controlling the dissociation dynam- 
ics by locally rebuilding appropriate phase space 
barriers. Furthermore, we hope that understand- 
ing the phase space dynamics of such higher di- 
mensional systems might provide a useful way of 
gaining insights into the nature of the fields com- 
ing out of an optimal control theory calculation. A 
recent work [73] addresses this issue in the context 
of a driven one degree of freedom system and the 
work presented in this paper is a first step towards 
generalizing to driven systems with higher degrees 
of freedom. 

5. One aspect of the problem that we have been silent 
about is that of the phenomenon of Arnold diffu- 
sion. In principle, the system considered here can 
undergo Arnold diffusion. However, in practice, the 
timescale for Arnold diffusion might be too long 
compared to the dissociation timescale and hence 
play no role in the mechanism of slow dissociat- 
ing trajectories. Nevertheless, from a fundamen- 
tal perspective, the current system provides an- 
other model, apart from the driven coupled quartic 
oscillator model[7S] considered by Chirikov in his 
seminal work, for studying Arnold diffusion and its 
quantum manifestations [ToT478] . 
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